Lorenz function of Bi2Te3/Sb2Te3 superlattices 
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Combining first principles density functional theory and semi-classical Boltzmann transport, the anisotropic 
Lorenz function was studied for thermoelectric Bi2Te3/Sb2Te3 superlattices and their bulk constituents. It was 
found that already for the bulk materials Bi2Te3 and Sb2Te3, the Lorenz function is not a pellucid function on 
charge earner concentration and temperature. For electron-doped Bi2Te3/Sb2Te3 superlattices large oscillatory 
deviations for the Lorenz function from the metallic limit were found even at high charge carrier concentrations. 
The latter can be referred to quantum well effects, which occur at distinct superlattice periods. 

PACS numbers: 31.15.A-,71.15.Mb,72.20.Pa,72.20.-i 
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INTRODUCTION 

For many decades thermoelectric (TE) energy conversion 
successfully enabled self-supporting energy devices for outer- 
space missions or integrated electronics QjJ, |2f] . However, bad 
conversion efficiency prohibited thermoelectric s the break- 
through as an alternative energy source. The conversion per- 
formance of a TE material is quantified by the figure of merit 
(FOM) 



ZT = 



-T = 



(1) 



where a is the electrical conductivity, S the thermopower, 
Kei = LoT and KW, are the electronic and lattice contribution 
to the thermal conductivity, respectively. L denotes the Lorenz 

function, which becomes the Lorenz number Ln = \? b ) in 
the highly degenerate, metallic limit. 

In recent years nano-structuring concepts ^ 01 enable 
higher values for ZT by increasing the numerator, called 
power factor PF = aS 2 , or decreasing the denominator of 
Eq. Q] The latter is obtained by phonon-blocking at super- 
lattice (SL) interfaces or grain boundaries fl5|-|8J and leads to 
a reduced lattice thermal conductivity K p i,. Here, the Lorenz 
function is particularly important for thermoelectrics, provid- 
ing a measure to separate the electronic and lattice contribu- 
tion to the thermal conductivity 19:]. Deviations L^= Lq already 
occur in the degenerate limit for simple metals, semi-metals 
and semi-conductors ifioll . Hence, assuming incorrect values 
for the Lorenz number leads to incorrect values for jq; and 
Kph and can even sum up to non-physically negative values 
for K p h llin . To the best of our knowledge, investigations on 
the Lorenz function of thermoelectric SLs on an ab initio level 
are missing so far. 

In the present work we analyse the anisotropic Lorenz func- 
tion for Bi2Te3/Sb2Te3 SLs, as well as for the bulk con- 
stituents. The two telluride single crystals and the composed 
p-type SL show highest values for bulk and nano-structured 
TE so far 01211 . On the basis of ab initio density functional 
theory (DFT) and semi-classical Boltzmann transport equa- 



tions (BTE) the Lorenz function is in particular studied for dif- 
ferent charge carrier concentrations and SL periods at room- 
temperature. 

METHODOLOGY 

For both Bi2Te3 and Sb2Te 3 , as well as for the com- 
posed SLs, we used the experimental lattice parameters and 
relaxed atomic positions [13] as provided for the hexago- 
nal Bi2Te3 crystal structure. The 15 atomic layers per unit 
cell are composed out of three quintuples Tei-Bi-Teo-Bi-Tei. 
The hexagonal lattice parameters are equally chosen to be 



^hex 
l BiTe 



4.384A and c%fr e 



30.487A for Bi 2 Te 3 , Sb 2 Te 3 and 
the SLs respectively. Preceding studies revealed that a larger 



in-plane lattice constant, e.g. a h %* Je > a ^L , is favourable for 



Jiex 



an enhanced cross-plane TE transport 11141 11511 . To introduce 
SLs with different SL periods we subsequently substitute the 
Bi sites by Sb, starting with six Bi sites in hexagonal bulk 
Bi 2 Te3. Substituting four atomic layers of Bi with Sb leads to 
a (Bi2Te3) x /(Sb 2 Te3)i_ v SL with x = 1, that is one quintuple 
Bi2Te3 and two quintuple Sb 2 Te3. The latter case coincides 
with a (10A/20A)-(Bi 2 Te3/Sb 2 Te 3 ) SL in the experimental no- 
tation ofRef. iQj] . 

Semi-classical BTE were extensi vely used in the past to cal- 
culate TE transport properties II16[ Il7ll and offer a high relia- 
bility for narrow-gap semi-conductors in a broad doping and 



temperature range 11141 118M20II . Within the relaxation time ap- 
proximation (RTA) the transport distribution function (TDF) 
2lll and with this the generalized conductance 



moments -Sf[j (M> ^) are defined as 



e=e: 



(2) 
de- 



f(fi t T) is the FERMl-DlRAC-distribution and v£^, 
note the group velocities in the directions in the hexago- 
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cross-plane directions are defined as 
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FIG. 1. Lorenz function L (thick black lines, ref. to left scale) 
and electronic contribution K e \ to the total thermal conductivity (thin 
green lines, ref. to the right scale) in dependence on position of 
the chemical potential fl within a spherical two band model. Re- 
sults are shown for (a) fixed effective masses m v j, = m c b and vary- 
ing temperatures and (b) fixed temperature T = 300 K and vary- 
ing effective masses. The band gap is fixed to E g = 0. 1 eV (gray 
shaded areas) and the Lorenz function is related to the metallic limit 
L = 2.44x 10" 8 Wfl/K 2 . 



nal basal plane and perpendicular to it, respectively. Within 
here the group velocities were obtained as derivatives along 
the lines of the Blochl mesh in the whole Brillouin zone 
(BZ) HI. The band structure Ej. of band v was ob- 
tained by accurate first principles density functional theory 
calculations (DFT), as implemented in the fully relativistic 
screened Korringa-Kohn-Rostoker Greens-function method 
(KKR) (22]. Within this approach the DlRAC-equation 
is solved self-consistently and with that spin-orbit-coupling 
(SOC) is included. Exchange and correlation effects were 
accounted for by the local density approximation (LDA) 
parametrized by Vosko, Wilk, and Nusair I123ll . Detailed stud- 
ies on the electronic structure, the thermoelectric transport and 
challenges in the numerical determination of the group veloc- 
ities of Bi?Te3, SbaTe.q and their SLs have been published be- 
fore (l4l[l5il24"" 



I25H . For convenience, the relaxation time T 
was chosen as lOfs for the considered systems. 

Straight forward, the temperature- and doping-dependent 
electrical conductivity a and thermopower S in the in- and 



<j ± ,,=e 2 J?l°l(Ll,T) 



1 3 



and the electronic part to the total thermal conductivity ac- 
counts to 



11 



i(^Mi(M,r) 



(4) 



The second term in eq. |4] introduces corrections due to the 
Peltier heat flow that can occur when bipolar conduction takes 
place i26l, l27ll . Using Eqs.|3]and|4]and the abbreviation K° = 



(5) 



^Jzf[ 2 ,j (jJ,, T) l2lll . we find the Lorenz function as 



Eq.|5]clearly shows that in the low temperature regime L con- 
sists of a constant term and a negative term of order T 2 . 



RESULTS 

To introduce our discussions, in Fig.Q~|the Lorenz function 
L and the corresponding electronic thermal conductivity K e i in 
dependence on the chemical potential jj. are shown for a spher- 
ical two band model (SBM). Varying temperatures and m c t, = 
m v }j (cf. Fig. Q2a)) and different effective masses ratio m cbjm vh 
and fixed temperature T = 300 K (cf. Fig.Ufb)) are assumed. 
Within here m c i, and m v b are the isotropic effective masses 
of the conduction band (CB) and valence band (VB), respec- 
tively. Setting the valence band maximum to zero and E g the 
band gap size, the TDF scales as «5f^(^,0) - % /^(-^i) 3 / 2 

and 3$}(n,0) - jm^,{\l - E g fl 2 for the VB and CB, re- 
spectively. From Eqs.|4]and|5]it is obvious that within a SBM 
deviations for L and K e i from the metallic limit will merely 
occur near the band gap, where the thermopower S changes 
significantly. Near the band edges S increases approximately 
as S ~ . Thus L, as well as K e \ minimize and the minimum 
decreases with decreasing temperature, while shifting towards 
the middle of the gap (cf. Fig.CQa)). At T = 100K - 0.8 
at the band edges. In the intrinsic regime and K e i in- 
crease, as the thermopower and electrical conductivity are 
reduced due to bipolar contributions. Figuratively speaking, 
the additional contribution arises from the fact that electron 
and holes can move together in the same direction, transport- 
ing energy but not carrying any net charge |]26|. According 
to Goldsmid 112811 and Price 12911 the deviation of the Lorenz 
number from the metallic limit in the intrinsic regime holds 
to some extent l/l q = 1 + j. m cb m vb (E g / kBT+A f, Therefore, 



2 {m cb +m rb y 

assuming a fixed charge carrier concentration, L / L o achieves 
very large values at small temperatures and/or large band 



3 



gaps. Assuming the above approaches 1281, 12911 . together with 



m c i, = m v b and E g = 0.1 eV one achieves l /l ~ 9 at room 
temperature for n located deep in the gap. If m v j, > m^, as 
shown in Fig.[TJb), the intrinsic regime N„ = N„ and with that 
the maximal value of l /lq and K e/ at bipolar conduction shifts 
towards the CBM. With increasing m v \, and hence due to the 
enhanced electrical conductivity a in the VB it is obvious, 
that K e \ under hole doping will increase, too. 

Fig. [2] presents first principle calculations for the Lorenz 
function L and the related electronic part K e \ of the thermal 
conductivity. The dependence on the charge carrier concen- 
tration for (a) Sb 2 Te 3 , (b) a (Bi 2 Te 3 ) x /(Sb 2 Te3)i- x SL at 
x = 2 /6 and (c) Bi 2 Te3, is shown, respectively. Due to the 
high conductivit y an isotropy a \\/a ± > 1 for all of the consid- 
ered svstems |[l4ll25ll . K e j j_ is strongly suppressed compared to 
K^H , too. Furthermore, it is obvious that the maximal peak of 
the Lorenz function is shifted towards the CBM, latter stem- 
ming from a larger density of states at the VBM and a higher 
absolute hole electrical conductivity. Maximal numbers for 
the Lorenz function i/f-o in the intrinsic regime were found 
to be between 6 and 10 for the considered systems, showing 
only a slight directional anisotropy. For Sb 2 Te3 (cf. Fig. [2 a)) 
the Lorenz function exhibits only minor anisotropics l \\/l ± in 
a wide doping range, while stating L± ~ 1.15L| at increased 
electron doping. Reduction of l h/l± due to bipolar diffusion 
effects is more apparent at hole doping compared to electron 
doping, here showing in-plane l /lq ~ 0.75 and l /l ~ 0.92 at 
an electron and hole doping of N = 3 x 10 19 cm~ 3 , respec- 
tively. For bulk Bi 2 Te3 the picture is comparable. However, 
in the thermoelectric most interesting range, about 200 me V 
around the band edges, l \\/l ± is always less than unity, com- 
parable to previous publications I20J. Furthermore Lj_ is 
larger than the metallic limit Lq. Very often values of L± ~ 
0.5 — 0.6 |L2, 3(J are assumed for the experimental determi- 
nation of K e / j_ in Bi 2 Te3/Sb 2 Te3 SLs. In turn, this most prob- 
ably leads to an underestimation of the electrical contribution 
to total thermal conductivity in cross-plane direction. For bulk 
Bi 2 Te3 experimental values 112611 for the in-plane part k^j are 
available as a reference in Fig. |2jc) (green, open circles). We 
find very good accordance to our calculations in the intrinsic 
range, while our results slightly overestimate K e n\ in the ex- 
trinsic regime. 

Strong deviations for the Lorenz function from the 
bulk limit could be found for an electron conducting 
(Bi 2 Te 3 y(Sb 2 Te 3 )i_ x SL at x = 2/e, i.e. (10A/20A)- 
(Bi 2 Te3/Sb 2 Te3). In a current publication ||25T| we showed, 
that strong quantum well effects (QWE) in the CB of the SLs 
lead to an enhanced electrical conductivity anisotropy, The 
latter was most pronounced for the SL at x = 2 /6 showing 
a \\/a ± ~ 20 at electron doping of N = 3 x 10 I9 cm~ 3 . Caused 
by the QWE, the cross-plane electrical conductivity a± is 
drastically suppressed, and hence L± remarkably enhanced. 
The cross-plane Lorenz function obtains rather large values 
between l /lq ~ 1 .5 — 1 .8 at extrinsic carrier concentrations of 
about N = 3 — 30 x 10 19 cm~ 3 for hole and electron doping, 
respectively. Additionally, oscillations of with varying 
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FIG. 2. Lorenz function L (solid lines, ref. to left scale) and elec- 
tronic contribution K e [ to the total thermal conductivity (dashed lines, 
ref. to right scale) in dependence on position of the chemical poten- 
tial for (a) bulk Sb 2 Te 3 (b) (Bi 2 Te 3 ) A -/(Sb 2 Te 3 ) 1 _ x SL at x= § 
and (c) bulk Bi 2 Te 3 . The in-plane (thick lines) and cross-plane (thin 
lines) transport directions are compared. The Lorenz function is re- 
lated to the metallic limit L = 2.44 x 1(T 8 WQ./K 2 . Plotted on to 
the graph of the Lorenz function in the in-plane direction is a color 
code referring to the charge carrier concentration. The red cross 
emphasizes the change from n to p doping. The temperature was 
fixed to 300 K. Thin vertical dash-dotted lines emphasize the posi- 
tion of the chemical potential for a charge carrier concentration of 
N = 3 x 10 19 cm~ 3 under p and n doping (red and blue color). The 
grey shaded areas show the band gap. Green open circles in (c) show 
experimental results from Ref. 02611 for K e [ || for an n-type Bi 2 Te 3 sin- 
gle crystal. 
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FIG. 3. Cross-plane component of the Lorenz function for 
(Bi2Te3) A -/(Sb2Te3)i_ A superlattices in dependence on the super- 
lattice period. The temperature is fixed to 300 K and results for 
three different charge carrier concentrations (in units of cm~ 3 ) are 
compared, (a) refers to electron doping, while (b) refers to hole 
doping. The Lorenz function is related to the metallic limit Lo = 
2.44 x 10~ 8 Wfi/K 2 . Lines are guides to the eye. 



doping are found, which are much more pronounced than in 
the bulk materials. Both effects have been proposed within a 
1 -dimensional model for thermoelectric SLs before Bill . As 
expected, in the extrinsic region, at increasing charge carrier 
concentration, L saturates gradually towards the metallic limit 
Lo and the thermal conductivity rises with electrical conduc- 
tivity. 



To support our findings, in Fig. [3] the cross-plane Lorenz 
function L± at different SL periods is shown. The influence 
of varying (a) electron and (b) hole doping is given, respec- 
tively. Under hole doping, due to the vanishing band-offset 
at the VBM |25[], the Lorenz function behaves as smooth in- 
terpolation between the bulk limits (cf. Fig. Ob)). At lower 
charge carrier concentrations l ±/lq is more suppressed due to 
a stronger impact of the bipolar diffusion. Under varying elec- 
tron doping we find l l/l being remarkably enhanced for SL 
periods of x — 2 /6 and x = 4 /6, respectively. For those SL peri- 
ods large suppressions of the cross-plane electrical conductiv- 
ity Oj_ were found, too. At N = 3 x 10 19 cm~ 3 anisotropics as 
large as a \\/a L ~ 20 and a \\/a ± ~ 14 for SL periods of x = 2 /6 
and x = 4 /6 are reported, respectively [|25ll ■ We state that due to 
quantum confinement effects in the electron-conducting SLs 
unexpected deviations of the Lorenz function from Lo can oc- 
cur also at higher values of doping, which is counterintuitive. 
The latter could lead to wrong estimations for the electronic 
part of the thermal conductivity K e / and consequently for the 
lattice thermal conductivity K p h- 



CONCLUSION 



We presented first principles calculation for the 
Lorenz function of electron- and hole-conducting 
Bi2Te3/Sb2Te3 superlattices and the related bulk materi- 
als at varying charge carrier concentration. As expected, due 
to bipolar conduction, the Lorenz function increases to large 
values within the intrinsic doping regime. More significantly, 
the Lorenz function L of the superlattices does not change 
monotonically at extrinsic charge carrier concentrations. 
While at increased doping an asymptote convergence of L 
towards the metallic limit Lo is found, a distinct oscillatory 
behaviour of L is observed. This is most pronounced under 
electron doping and caused by quantum well effects in the 
conduction bands of the superlattices. This counterintuitive 
effect has consequences for the determination of the thermal 
conductivity, as L is generally used to separate JQ/ and K p ] v 
At thermoelectrically profitable charge carrier concentrations 
the application of the metallic value Lq to determine the 
electronic thermal conductivity could lead to a deviation of 
a factor of two in both directions in the worst case. Conse- 
quently, this leads to wrong estimations of the lattice thermal 
contribution and the figure of merit. A similar behaviour was 
found theoretically for p-type SiGe superlattices Il32ll and 
this behaviour could be a general effect in thermoelectric 
superlattices influenced by quantum well effects. 
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